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Probing the dark side: Constraints on the dark energy equation of state from CMB, 

large scale structure and Type la supernovae 
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We have reanalysed constraints on the equation of state parameter, wq = P/ p, of the dark 
energy, using several cosmological data sets and relaxing the usual constraint vjq > — 1. We find 
that combining Cosmic Microwave Background, large scale structure and Type la supernova data 
yields a non-trivial lower bound on Wq. At 95.4% confidence we find, assuming a flat geometry of 
the universe, a bound of —2.68 < wq < —0.78 if wq is taken to be a completely free parameter. 
Reassuringly we also find that the constraint wq > — 1 does not significantly bias the overall allowed 
region for wq. When this constraint is imposed the 95.4% confidence bound is — f < wq < —0.71. 
Also, a pure cosmological constant (w = —1) is an excellent fit to all available data. Based on 
simulations of future data from the Planck CMB experiment and the SNAP and SNfactory supernova 
experiments we estimate that it will be possible to constrain wq at the 5 % level in the future. 

PACS numbers: 98.80.-k,12.60.-i,14.80.-j 



I. INTRODUCTION 

Several independent methods of observation all sug- 
gest that most of the energy density in the universe is in 
the form of a component with negative pressure - dark 
energy. The simplest possibility for such a component is 
a cosmological constant which has a constant equation 
of state Pq = —pq. However, in the general case the 
dark energy can have an equation of state which is time- 
dependent, Pq = w Q (t)p Q @ |,| g |,|, @ ||. Dark 
energy with a time-dependent equation of state has been 
invoked to explain the coincidence problem, the fact that 
the energy density in dark energy is roughly equal to that 
in dark matter exactly at the present epoch. By coupling 
a scalar field to matter one can obtain tracking solutions 
for the time dependence of the dark energy density so 
that it always follows the dominant energy density com- 
ponent. 

Generic to most of these proposed candidates for dark 
energy is the fact that wq > — 1 at all times. This is, e.g., 
the case for most quintessence models where a scalar field 
is rolling in a potential (potential driven quintessence). 
Since most of the plausible models lie in this category, 
the likelihood analyses which have been used to find wq 
have cut away the region with wq < — 1. From a purely 
phenomenological point of view this is not justified and 
can lead to severe bias in parameter determination. This 
is particularly worrisome because the most recent data 
set an upper limit on wq < —0.85 (68% confidence) 
(see also (jTl]), but not a lower limit. Therefore, if the 
true model has wq < —1, the upper bound could be 



'Electronic address: 
t Electronic address: 



tccn@nordita.dk 



:dvard@physto.sc 



wrong (and in principle a full analysis could even rule 
out a cosmological constant as being the dark energy). In 
fact there are several models which predict a dark energy 
component with w Q < -1 @ [H§ 0, ||, || This 
possible bias problem was also discussed in Ref. |18] . 

In the present paper we reanalyse cosmological data 
from the Cosmic Microwave Background (CMB), large 
scale structure (LSS) and Type la supernovae (SNe) 
without the constraint wq > — 1. We make the simpli- 
fying approximation that wq(€) = const. Even though 
this is certainly not true for many models of dark energy, 
almost all models can be very well approximated by a 
model having a constant WQ iC ff , which is then calculated 
as a properly weighted mean of wq (t) |l9| . With this 
approximation the dark energy density evolves simply as 
Pq cx a~ 3 ( 1+WQ \ where a is the scale factor. 

Our analysis of the present data is in many ways similar 
to that performed in Ref. |IC[ ]. However, in addition to 
the extension of the parameter space to wq < — 1 we also 
include new data from the 2dF galaxy survey. 

Finally, we discuss the prospects for measuring wq pre- 
cisely with future high precision CMB and SN data. 



II. DATA ANALYSIS 

A. CMB and large scale structure 

CMB data set — Several data sets of high precision 
are now publicly available. In addition to the COBE |2C[| 
data for small I there are data from BOOMERANG [glj, 
MAXIMA H|, DASI H and several other experiments 
p4] , p5[ . Wang, Tegmark and Zaldarriaga |24j have com- 
piled a combined data set from all these available data, 
including calibration errors. In the present work we use 
this compiled data set, which is both easy to use and 
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includes all relevant present information. 

LSS data set — At present, by far the largest survey 
available is the 2dF |^6| of which about 147,000 galaxies 
have so far been analysed. Tegmark, Hamilton and Xu 
p7| have calculated a power spectrum, P(fc), from this 
data, which we use in the present work. The 2dF data 
extends to very small scales where there are large effects 
of non-linearity. Since we only calculate linear power 
spectra, we use (in accordance with standard procedure) 
only data on scales larger than k — 0.2h Mpc -1 , where 
effects of non-linearity should be minimal. 

The CMB fluctuations are usually described in terms 
of the power spectrum, which is again expressed in terms 
of Ci coefficients as 1(1 + l)Ci, where 



Ci = (|a im | 2 ). 



(1) 



The ai m coefficients are given in terms of the actual tem- 
perature fluctuations as 



I in 



(2) 



Given a set of experimental measurements, the likelihood 
function is 



£(B) x < x]> | -Ix^dey^x 



(3) 



where 9 = (il, fit, H ,n,r, . . .) is a vector describing the 
given point in parameter space, a; is a vector containing 
all the data points, and C(6) is the data covariance ma- 
trix. This applies when the errors are Gaussian. If we 
also assume that the errors are uncorrelated, it can be 
reduced to the simple expression, C oc e~ x / 2 , where 
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is a x 2 -statistic and iV max is the number of power spec- 
trum data points (2{| . In the present paper we use Eq. (H) 
for calculating \ 2 f° r tne CMB data. Since we also use 
data from the 2dF survey the total % 2 is then given by 
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(5) 



The procedure is then to calculate the likelihood func- 
tion over the space of cosmological parameters. For cal- 
culating CMB and matter power spectra we have used 
the publicly available CMBFAST package ||. The 2D 
likelihood function for (£l m ,WQ) is obtained by keeping 
{0, m , wq) fixed and marginalising over all other parame- 
ters. 



As free parameters in the likelihood analysis we use 
fl m , the matter density, wq, the dark energy equation 
of state, fib, the baryon density, Hq, the Hubble pa- 
rameter, n, the scalar spectral index, and r, the optical 
depth to reionization. The normalisation of the CMB 
data, Q, and of the 2dF data, b, are taken as completely 
free and uncorrelated parameters in the analysis. This 
is very conservative and eliminates any possible system- 
atics involved in determining the bias parameter. We 
constrain the analysis to flat (f2 m + £lQ = 1) models, and 
we assume that the tensor mode contribution is negligi- 
ble {T/S = 0). These assumptions are compatible with 
analyses of the present data [Q , and relaxing them does 
not have a big effect on the final results. 



TABLE I: The different priors used in the analysis. 



Parameter 


Prior 




0.1-1 


WQ 


-3 - -0.5 


Q, b h 2 


0.020 ± 0.002 (Gaussian) 


h 


0.72 ± 0.08 (Gaussian) 


n 


0.5-1.4 


r 


0-1 


Q 


free 


b 


free 



Table 1 shows the different priors used. We use 
the constraint Hq — 72 ± 8kms _1 Mpc _1 [h = 
H /(W0 kms-iMpc" 1 )] from the HST Hubble key 
project Q (the constraint is added assuming a Gaus- 
sian distribution) and the constraint H,t,h 2 = 0.020±0.002 
from BBN ||. 

Fig. 1 then shows the 68.3% and 95.4% confidence 
allowed regions, corresponding to A% 2 = 2.31 and 6.17 
respectively. 

For very negative values of wq the CMB and LSS con- 
straint is independent of wq . The reason is that at such 
low wq the dark energy influences the CMB spectrum 
only via the late integrated Sachs- Wolfe effect. However, 
the late ISW effect decreases in magnitude as wq de- 
creases. This effect can be seen in Fig. 2 where CMB 
power spectra are plotted for different values of wq. 



B. Type la supernovae 

SN data set — The SN data set used in this analy- 
sis corresponds to Fit C from the Supernova Cosmology 
Project as described in j32|. This is a subsample of a 
total of 60 SNe where two SNe are excluded as statisti- 
cal outliers, two because of atypical lightcurves and two 
because of suspected reddening. 

We use this data set to fit fl m and wq, taking advan- 
tage of the cosmology dependence of the distance-redshift 
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FIG. 1: The 68.3 % (dark shaded) and 95.4 % (light shaded) 
confidence allowed regions for Q m and wq using CMB, HST, 
BBN and LSS data. 
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FIG. 2: Different CMB power spectra for different values of 
wq. In all cases the model parameters are those of the fiducial 
ACDM model, ft m = 0.3, Q A = 0.7, n b h 2 = 0.020, H = 
70 kms - Mpc -1 . The full line is for wq = —1, the dashed for 
wq — —2, and the dotted for wq — —4 (in order of decreasing 
Ci at low I). 

relation. Type la SNe are very useful as distance indi- 
cators because of their high luminosities and small dis- 
persion among their peak absolute magnitudes [a m ~ 
0.15). Also, they have distinct spectral lines, allowing 
for accurate redshift determinations. 

The apparent and absolute magnitudes are related by 

m{z) = M + 5 log [D L (z, Cl m , w Q )} - 5 logff + 25, (6) 

where T>l := HqcIl is the part of the luminosity distance 



that remains after multiplying out the dependence on 
the Hubble constant (expressed here in units of km s _1 
Mpc -1 ). In the low redshift limit, Eq. ([j|) reduces to a 
linear Hubble relation between m and log z : 



m(z) = M + 51ogz, 



(7) 



where we have expressed the intercept of the Hubble line 
as A4 := M — 51ogiJ + 25. This quantity can be mea- 
sured from the apparent magnitude of low redshift stan- 
dard candles, without knowing the value of Hq. Thus, 
with a set of apparent magnitude and redshift measure- 
ments m(z) for Type la SNe, we can find the best-fit 
values of (Q m ,WQ) to solve the equation 

m(z)-M = 5log[D L (z,n m ,w Q )}. (8) 

The x 2 is then given by 



E 



5 log [V L (z t ,n m ,w Q )}-MY 



(9) 



where Oi is the statistical uncertainty for each event, We 
assume a flat geometry of the universe when calculating 
T>l and marginalise over M. where we assume no prior 
knowledge. 

In Fig. |, the 68.3% and 95.4% confidence allowed 
regions are showed. The best-fit values are ft m = 0.45 
and wq = —1.9, indicating the possibility of a bias in the 
parameter determination when imposing the constraint 
wq > -1. 
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FIG. 3: The 68.3% (dark shaded) and 95.4% (light shaded) 
confidence allowed regions for fi m and wq using the 54 type 
la SNe from the Supernova Cosmology Project. 



C. Combined constraint 

When all the available data is combined, a fairly strin- 
gent bound on wq is obtained. The 68.3% and 95.4% 
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confidence combined bounds are shown in Fig. 4. For 
wq alone we find a 95.4% confidence bound of —2.68 < 
w Q < -0.78. 
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FIG. 5: The 68.3% (dark shaded) and 95.4% (light shaded) 
confidence allowed regions for £l m and wq using all available 
data and imposing the bound wq > — 1. 



FIG. 4: The 68.3 % (dark shaded) and 95.4 % (light shaded) 
confidence allowed regions for £l m and wq using all available 
data. 



III. DISCUSSION 
A. Constraints from present data 

An important point is to determine whether the re- 
laxation of the bound wq > — 1 significantly affects the 
likelihood analysis for the part of parameter space which 
is above wq = — 1. This could be the case if the best-fit 
value of wq lies in the excluded region as is the case for 
the current SN data, where the best-fit value corresponds 
to Q m = 0.45, wq = -1.9. For the CMB+LSS data the 
situation is the same, with the overall best fit being at 
n m = 0.26, wq — —2.6. In Fig. 5 we plot the the like- 
lihood contours with the constraint wq > — 1 imposed. 
In terms of wq alone the 95.4% confidence bound is now 
-1 < w Q < -0.71. 

Comparing this to the bound obtained without the 
constraint shows that the likelihood contours are fairly 
similar and that the bias is not a significant problem. 
Our Fig. 4 can for instance also be compared to Fig. 3 of 
Rcf. 0. Apart fr om the smaller allowed region in O m 
because of the tighter BBN constraint and the use of the 
full 2dF data set, the two plots are very similar. The 



bound on wq found in Ref. 



It 



is 



-1 < w < —0.73, very 



similar to the constraint found in the present analysis 
when the constraint wq > — 1 is imposed. 

This shows both that our constrained analysis is con- 
sistent with Ref. jlO| and that this previous analysis of 
wq is not seriously biased. 



The next very important point of our analysis is that 
one obtains a non-trivial lower bound on wq from the 
combination of CMB, LSS and SN data. From SN 
data alone we can infer that wq > —12 is ruled out 
at the 68.3% confidence level, but combining this with 
CMB and LSS data tightens the bound significantly to 
wq > —2.68 at 95.4% confidence. It is also very interest- 
ing and perhaps somewhat suggestive that a cosmological 
constant lies in the 68.3% confidence allowed region. 



B. Constraints from future data 

The ability to constrain the equation of state param- 
eter wq of a dark energy component using future CMB 
and Type la SN data have been recently investigated by 
a number of authors, see, e.g., |33| and references therein. 
Our analysis differ in the respect that we do not impose 
the constraint wq > — 1 and that we use the most current 
anticipated data sets. 

CMB data set — For CMB we use simulated data from 
the Planck Surveyor satellite. For simplicity we use 
only data from the HFI 100 GHz channel, assuming an 
angular resolution of 10.7 arcmin and a pixel noise of 
AT/T = 1.7 x 10" 6 @. This channel is not polariza- 
tion sensitive and so our assumed data set seems con- 
servative compared to what can be expected from the 
full Planck data. On the other hand we do not include 
foregrounds in our analysis. The simulated data is gen- 
erated from an underlying flat model with the follow- 
ing parameters: Q m = 0.3, wq = —1, ili,h 2 = 0.02, 



H = 70 kms^Mpc" 



1.0, and r = 0. 



SN data set — We use simulated data sets correspond- 
ing to three year's data from the proposed satellite tele- 



5 



scope the Supernova/ Acceleration Probe (SNAP; |35 ) 
and the predicted results from the Supernova Factory 
Campaign (SNfactory; scheduled to start in the end 
of 2002. 

The SNAP satellite would be capable of discovering 
and taking spectra of ~ 2800 Type SNe per year for red- 
shifts z < 1.7. The current projected redshift distribu- 
tion follows the distribution in the SNAP proposal |3f| for 
z < 1 and is approximately uniform at higher redshifts 
jj7|- The SNfactory data set consists of 200 SNe between 
0.03 < z < 0.06 and 100 SNe between 0.06 < z < 0.15. 
The simulated data is generated for the same underlying 
model as for the CMB simulations. 
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FIG. 6: The 68.3 % (dark shaded) and 95.4 % (light shaded) 
confidence allowed regions for Q m and icq using simulated 
data from the Planck satellite and the SNAP and SNfactory 
observatories. 



roughly 0.05 relative precision in the two parameters. In 
this note we have neglected the possible use of multiple 
imaged core-collapse SNe to constrain wq and Q m , see 
Ref. ||. 

It is important to note that the confidence regions pre- 
sented only take into account statistical uncertainties. 
For future Type la SN data sets, systematic errors from, 
e.g., dust obscuration, luminosity evolution and gravi- 
tational lensing might be comparable to or even larger 
than the statistical errors. In Fig. 7 the confidence re- 
gions corresponding to the middle panel of Fig. 6 is shown 
if lensing effects from 90 % NFW dark matter halos and 
10 % point masses are included in the simulated data set 
(generated with Q m = 0.3 and wq = —1). It is obvious 
that gravitational lensing, if not taken into account, will 
cause an underestimation of the matter density and in 
order to reach the full potential of the large statistics, we 
need to to correct for the effect, e.g., as described in [ p9[ . 
For the CMB data there may also be systematic errors of 
a magnitude comparable to the purely statistical ones. 

However, if the systematic errors can be controlled in 
an effective matter, we conclude that it should be possible 
to constrain the equation of state parameter wq of the 
dark energy to high accuracy using future CMB and SN 
data sets. 
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FIG. 7: The 68.3% (dark shaded) and 95.4% (light shaded) 
confidence allowed regions for Q m and wq using simulated 
data from the SNAP and SNfactory observatories, with lens- 
ing effects from 90% NFW halos and 10% point masses in- 
cluded. The simulated data is generated with Q m = 0.3 and 

WQ = -1. 



In Fig. 6 we show a likelihood analysis based on the 
simulated CMB and Type la SN data. It is clear that 
again the data will complement each other, the CMB 
data being very sensitive to O m and the SNIa data to 
wq. From the combined data set we estimate that it is 
possible to obtain a 95.4% confidence interval on wq of 
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